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ABSTRACT 

We present here the first results from the Chandra ERA (Environments of Radio-loud AGN) Large Project, 
characterizing the cluster environments of a sample of 26 radio-loud AGN at 2 ~ 0.5 that covers three decades 
of radio luminosity. This is the first systematic X-ray environmental study at a single epoch, and has allowed 
us to examine the relationship between radio luminosity and cluster environment without the problems of 
Malmquist bias. We have found a weak correlation between radio luminosity and host cluster X-ray luminos- 
ity, as well as tentative evidence that this correlation is driven by the subpopulation of low-excitation radio 
galaxies, with high-excitation radio galaxies showing no significant correlation. The considerable scatter in the 
environments may be indicative of complex relationships not currently included in feedback models. 
Subject headings: galaxies: active, galaxies: jets 



1. INTRODUCTION 

Understanding how the properties of radio-loud AGN relate 
to their cluster environments is crucial for our understanding 
of the role of AGN feedback in galaxy evolution. Suppres- 
sion of star formation by feedback from this type of AGN is 
now an important feature of simulations of galaxy evolution, 
allowing them to match the observed galaxy sizes and star for- 
mation rates (Croton et al. 2006). The fact that the radio jets 
disturb the cluster environment can be seen in detailed stud- 
ies of nearby radio galaxies (eg Kraft et al. 2003; Fabian et al. 
2003; Forman et al. 2005; Croston et al. 2009). The processes 
involved are described in the review by McNamara & Nulsen 
(2007). Galaxy feedback via radio jets is however a complex, 
two-way process since the environment is also expected to af- 
fect the properties of the radio galaxy, and the relationship 
between energy input from the radio-loud AGN and environ- 
ment, and how their relationship evolves with epoch is as yet 
poorly understood. The radio jets transport energy a consid- 
erable distance into the cluster and are themselves modified 
by the intra-cluster medium. Do the properties of the large- 
scale cluster environment in their turn affect the feedback loop 
maintaining the AGN, or are the AGN properties determined 
by the more local environment of the host galaxy? And how 
does does this disruption of the cluster environment affect its 
evolution? 

Thus two outstanding questions are whether the radio lu- 
minosity is related to the large-scale cluster environment, and 
whether typical environments evolve with epoch. 

These questions have been in consideration for some time. 
At low redshifts, it has long been known that FRI galaxies 
appear to inhabit richer clusters than the more luminous FRII 
galaxies (eg Longair & Seldner 1979; Prestage & Peacock 
1988). However, Prestage & Peacock found that although 
there was a clear difference in average richness, there was also 
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a large scatter in the richness of environments, with both types 
having examples at similar extremes. Hill & Lilly (1991) ex- 
tended this work by comparing the low redshift results with a 
flux-limited sample at z ~ 0.5. They found that at the higher 
redshift range, the FRII galaxies were spread more evenly 
over a wider range of cluster richnesses than at low redshift, 
raising the possibility that the environments of FRII galaxies 
evolve. 

There are a number of selection biases that may have af- 
fected these early studies, and that still make sample selection 
difficult today. The overwhelming majority of radio galaxies 
in the local universe are low luminosity FRI galaxies. How- 
ever, with flux-limited samples we see only increasingly lu- 
minous objects, predominantly FRIIs, as we increase redshift 
- Figure [T| (left) shows the problem clearly. This gives op- 
posite biases in results at low and high redshifts and makes 
systematic comparisons between low and high redshift data 
difficult. 

There is also a potential confounding factor in the choice 
of subsamples for comparisons. The studies cited above use 
the FR classification (Fanaroff & Riley 1974), which is based 
on the location of the brightest area of luminosity in the radio 
lobes. Fanaroff & Riley found that this corresponded to a rea- 
sonably clear division in radio luminosity, and this division 
was later found also to be related to optical luminosity (Owen 
& Ledlow 1994). However, radio-loud AGN can also be clas- 
sified by their optical spectral properties, ie whether or not 
the object has strong emission lines (eg Hine & Longair 1979; 
Laing et al. 1994). High-excitation radio galaxies (HERGs), 
which have strong emission lines, incorporate both broad- and 
narrow-line radio galaxies (BLRGs and NLRGs), while low- 
excitation radio galaxies (LERGs) lack strong emission lines. 
The majority of HERGs are luminous FRII galaxies, whereas 
LERGs span the full range of radio luminosities, including al- 
most all FRIs and a significant number of FRIIs. Hine & Lon- 
gair, for example, found that about 10% of weak radio galax- 
ies were HERGs, with the percentage increasing to over 70% 
at high radio luminosities. Willott et al. (2001) also found 
that the proportion of HERGs increased with radio luminos- 
ity, with examples of HERGs and LERGs at all luminosities. 
Thus subpopulations based on the Fanaroff-Riley classifica- 
tion will consist of mainly, but not exclusively, LERGs for the 
FRI class and a mixture of HERGs and LERGs for the FRII 
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class. 

Discussions of AGN unification (eg Antonucci 1993; Laing 
et al. 1994) recognized that although radio-loud QSOs, 
BLRGs and NLRGs could be interpreted as the same class 
of object viewed at different orientations, LERGs appeared to 
be a different class of object. Further differences between the 
two classes emerged (eg Whysong & Antonucci 2004; Ogle 
et al. 2006; Chiaberge et al. 2002; Evans et al. 2006; Hard- 
castle et al. 2006) and it is now thought that HERGs are un- 
dergoing radiatively efficient accretion, while LERGs are ra- 
diatively inefficient (eg Hardcastle et al. 2007). The classes 
therefore appear to be fundamentally different. 

More recent studies into the effects of cluster richness there- 
fore tend to split samples by spectral class rather than, or as 
well as, FR class. Looking at relatively low redshift samples 
(z < 0.4), Best (2004) and Hardcastle (2004) both found that 
the different types of HERG all inhabit poor environments, 
but that LERGs are spread across a broader range of environ- 
ments. Furthermore, the environments of the different types 
of HERG were all similar, as would be expected if they were 
different aspects of the same class of object. Gendre et al. 
(2013) obtained the same results for HERGs and LERGs for 
z < 0.3, and also found that the FR class was independent of 
excitation class. At higher redshifts, Harvanek et al. (2001) 
found QSOs at 0.4 < z < 0.65 in richer environments than 
at z < 0.4, and Belsole et al. (2007), using a sample of FRII 
galaxies that were mostly HERGs (0.45 < z < 1.0), also 
found them inhabiting relatively rich environments. Since 
HERGs are mostly FRII galaxies and LERGs are both FRI 
and FRIIs, these results are compatible with Hill & Lilly's 
results for FRI and FRII galaxies. 

Does this imply that the environments of HERGs change 
with time, or, given the combination of the Malmquist bias 
and the paucity of high luminosity radio galaxies at low red- 
shifts, does this imply that the environments of high luminos- 
ity sources are typically richer than those of low luminosity 
sources? Best (2004), using a sample with z < 0.1, found 
a strong correlation between radio luminosity and environ- 
ment richness for LERGs, but not for HERGs. Belsole et al. 
(2007) found no correlation for their high-z FRII sample, but 
Wold et al. (2000; 0.5 < z < 0.9) and Falder et al. (2010; 
z ~ 1.0) both found a correlation for high redshift, radio-loud 
QSOs. There is therefore some evidence that radio luminos- 
ity is related to the richness of the cluster environment for 
at least some classes of radio galaxy, but this does not ex- 
clude the possibility of evolution with epoch. However, Wold 
et al. (2000) and McLure & Dunlop (2001) both compared 
compared their samples with results from studies at different 
redshifts and found no evidence of a variation of environment 
with redshift - evidence supporting a link between environ- 
ment and radio luminosity rather than epoch. So the picture 
of the relationship between radio-loud AGN and their large- 
scale cluster environments, and for the evolution of those en- 
vironments, remains confused, and the long-term aim of our 
research program is to clarify these issues. 

The studies described above use a variety of measures of 
cluster richness, raising questions about the compatibility of 
their results, and how well the total cluster mass is being 
traced. The scaling relations between mass, luminosity and 
temperature of the ICM are well defined for galaxy clusters 
(eg Arnaud et al. 2007; Vikhlinen et al. 2006; Pratt et al. 
2009), and so the X-ray luminosity of the ICM is one of the 
most well understood measures of cluster richness. Despite 
the scatter in the mass-luminosity relationship, it should pro- 



vide a more accurate proxy for cluster mass than the optical 
measures. There is evidence that it is strongly related to the 
optical measures; for example, Yee & Ellingson (2003) have 
found a correlation between the galaxy cluster center corre- 
lation amplitude B gc and X-ray luminosity, and Ledlow et 
al. (2003) between Abell class and X-ray luminosity. X- 
ray observations also have the advantage of directly probing 
the medium into which the radio jets are propagating, and so 
the relationship between X-ray cluster properties and radio- 
source characteristics might be expected to be tighter than for 
measures associated with the cluster galaxy population. 

In this paper, we present the first results of the Chandra 
ERA (Environments of Radio-loud AGN) program, whose 
aim is to make a systematic, X-ray based examination of the 
effects of cluster environment and epoch on the properties 
of radio-loud AGN for the first time. If the relationship be- 
tween these properties can be clearly disentangled, then well- 
characterized relationships could be used for modeling radio- 
jet feedback in galaxy evolution models. In order to separate 
the effect of luminosity from that of redshift, we use a sample 
limited to a narrow redshift range but spanning a wide range 
of radio luminosities. This paper contains the results of this 
first part of this program. In future papers, we will compare 
these results with observations at different epochs to look for 
evidence of environmental evolution with redshift. 

Throughout this paper we use a cosmology in which Hq = 
70 km s _1 Mpc~\ tt m = 0.3 and S1 A = 0.7. Unless other- 
wise stated, errors are quoted at the ler level. 

2. THE SAMPLE 

We made use of the sample of McLure et al. (2004), 
which contains all narrow-line and low-excitation radio galax- 
ies with redshifts between 0.4 and 0.6 from four flux-limited, 
spectroscopically complete, low-frequency radio surveys of 
the northern hemisphere - 3CRR (Laing et al. 1983), 6CE 
(Eales et al. 1997; Rawlings et al. 2001), 7CRS (Lacy et 
al. 1999; Willott et al. 2003), TexOx-1000 (Hill & Rawl- 
ings 2003). This sample is ideal for achieving our aim of 
comparing radio luminosity and environment richness with- 
out contamination by evolution, as it covers three decades in 
radio luminosity in a small redshift range while being dis- 
tant enough to contain high luminosity radio galaxies but near 
enough for low luminosity galaxies still to be detectable and 
for X-ray observations of the ICM to be feasible with reason- 
able exposure times. 

The luminosity and redshift range for the four surveys are 
shown on the left of Figure Q] with the McLure et al. (2004) 
sample highlighted. The sample contains 41 sources cover- 
ing three decades of radio luminosity, and includes high- and 
low-excitation sources and a range of radio source morpholo- 
gies. Host galaxy properties for the full sample were derived 
by McLure et al. (2004) from HST WFPC2 observations, 
and an optical environment study has recent been carried out 
by Herbert et al. (2013), using the spatial covariance func- 
tion B gq . Because of the amount of X-ray observing time 
required, we could not use the entire McLure et al. sample 
for this project. We therefore constructed the ERA sample, 
a representative subsample of 24 sources from the McLure et 
al. sample covering its full range of radio luminosity and con- 
taining the same subgroups. We limited our sample to sources 
with extended lobes, so compact sources were excluded from 
our selection. Ten of the sources had already been imaged 
in X-ray by XMM-Newton or Chandra, and we obtained new 
Chandra and XMM-Newton observations of fourteen sources 
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for this project. Subsequently, X-ray data for another two 
sources from the parent sample have become available and 
they have been included. Figure Q] shows radio luminosity 
plotted against redshift for the ERA sample (right) alongside 
the McLure et al. (2004) sample (left), and Figure |2] shows 
the coverage of the ERA sample. 

Table Q] lists details of the 26 sources. Positions and red- 
shifts were taken from McLure et al. (2004). For all galax- 
ies except 3C 457, we obtained Galactic column densities 
from Dickey & Lockman (1990) via the HEASARC tools; for 
3C 457, we used the higher column density found by Konar 
et al. (2009). Excitation type was taken from McLure et al. 
(2004), except for 3C 295. 3C 295 is classified as a LERG in 
the on-line 3CRR catalog based on the results of Lawrence et 
al. (1996), but this classification is questioned (Varano et al. 
2004). We have here classified 3C 295 as a LERG, but where 
appropriate have analyzed LERG subsamples with and with- 
out 3C 295. 7C 021 9+3423 's classification is also uncertain; 
we have followed McLure et al. in classifying it as a possible 
HERG. 

3. OBSERVATIONS AND DATA PREPARATION 
3.1. X-ray data 

The X-ray observations for this study came from Chandra 
and XMM-Newton. The 3C sources had already been ob- 
served; the 6C, 7C and TOOT observations were made for 
this program. The XMM-Newton observations used the three 
EPIC cameras with the medium filter and the Chandra obser- 
vations used the ACIS-S3 chip in either FAINT or VFAINT 
mode. Observation IDs and times are given in Table |2] 

We used the Chandra analysis package CIAO v4.3 for pro- 
cessing the Chandra events file. We reprocessed the files us- 
ing the chandrajrepro script, applying particle background 
cleaning for observations in the VFAINT mode, and then 
checked for background flares by extracting a light curve us- 
ing dmextract. We excluded events at more than 3er above 
the background rate - we removed a large flare from 3C 295 
and small flares from 3C 16, 3C 200, 7C 0219+3423 and 
TOOT 1626+4523. Screened observation times are included 
in Table |2] 

The XMM-Newton events files were reprocessed with the 
latest calibration data using XMM-Newton SAS vl 1.0.0. The 
pn camera data were filtered to include only single and dou- 
ble events (PATTERN < 4) and data from the MOS cameras 
were filtered using the standard pattern mask (PATTERN < 
12). The data sets were also filtered to remove bad pixels, 
bad columns etc. We checked each events file for flares using 
the light curves at higher energy levels than those emitted by 
the sources (10-12 keV for the MOS cameras, 12-14 keV for 
the pn camera) and used good-time-interval (GTI) filtering to 
select data where the light curve was within 20% of the quies- 
cent level. We then used evigweight to correct the events files 
for vignetting. 

The particle background in the XMM-Newton sources 
was removed using the method described by Croston et al. 
(2008a). This uses closed filter files (supplied courtesy of 
E. Pointecouteau) that were processed, filtered and weighted 
in the same manner as the source data sets. The closed fil- 
ter data were rotated to match the source observations and 
scaling factors were calculated by comparing the count rates 
at 12-14 keV (pn camera) and 10-12 keV (MOS cameras). 



The closed filter data were then scaled by these factors before 
carrying out background subtraction when generating profiles 
and spectra. 

For both the Chandra and XMM-Newton sources, we ex- 
amined images of the sources in ds9, overlaying them with 
the radio contours to find the approximate center of the ex- 
tended source. We identified point sources in the data sets and 
emission associated with the radio lobes, which we excluded 
during subsequent analysis. 

3.2. Radio data 

Radio maps were used to make the overlay images shown 
in Appendix [A] and to mask out the radio lobes so that any 
radio-related X-ray emission did not contaminate our mea- 
surements of the cluster properties. In many cases we used 
existing maps, either from the 3CRR Atlasj; the Faint Im- 
ages of the Radio Sky at Twenty Centimeters Survey (FIRST) 
(Becker et al. 1995); Hardcastle et al. (2002b); and Cros- 
ton et al. (2005b). For some of the less luminous sources 
where adequate maps were not available from these sources, 
we used the 1.4-GHz observations of Mitchell (2005) which 
were taken in A and C configurations. We obtained these from 
the VLA archive (Program ID AR477) and reduced them in 
AIPS in the standard manner. Table [3] contains full details of 
the radio maps used. 

4. ANALYSIS 

The aim of the analysis was to find the temperature and 
X-ray luminosity of the ICM emission of the radio galax- 
ies. Where possible, the temperature was obtained by spec- 
tral analysis; when there were insufficient counts, it was esti- 
mated from the count rate. The luminosity was determined by 
integrating the su rface brightness profiles to the -R500 radius 
(defined in Section l4~2b . 

4.1. Imaging 

AppendixlAlcontains images of the X-ray emission of each 
cluster in the 0.5-5 keV energy range, with the radio emis- 
sion overlaid as contours. We generated images for the XMM- 
Newton sources using the method described in Croston et al. 
(2008a). An image was extracted for each of the three EPIC 
cameras using evselect. The MOS images were then scaled 
to make their sensitivity equivalent to the pn camera image so 
that there would be no chip-gap artifacts when the three im- 
ages were combined. We generated exposure maps for each 
camera using eexpmap, which were used to correct for the 
chip gaps, but not for vignetting as this leads to incorrect scal- 
ing of the particle background that dominates at large radii. 
The resulting images are therefore not vignetting corrected; 
they are purely pictorial and not used in any subsequent anal- 
ysis. 

For both Chandra and XMM-Newton sources, we used dm- 
filth to replace point sources by Poission noise at the level of 
nearby regions, and then applied Gaussian smoothing using 
aconvolve. We then used ds9 to display the X-ray emission 
and overlay the radio emission contours. 

4.2. Spatial analysis 

We extracted a radial surface brightness profile from the 
events file of each source by taking the average counts in 
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TABLE 1 

The ERA sample, in order of radio luminosity 



Source 


RA (J2000) Dec 
h m s ° ' " 


Redshift 


Scale 
kpc arcsec -1 


logio L151 
WHz" 1 sr- 1 


Radio 
structure 11 


Type 


nH 

xlO 20 cm- 2 


lOOl 1 .301 +.3658 


13 Ul 25.U3 


i i£. co no a 
+36 58 (jy.4 


0.424 


5.56 


24.76 


J 


LERG 


1.32 


1UU1 1255+3556 


n ?c cc ci 
12 jj jj.53 


. or r/r tc q 
+3j jo 3j.o 


f\ Alt 

U.4 / 1 


c m 
5.91 


oc m 
25.01 


J 


LERG 


1.30 


1UU1 1626+4523 


1 f. If, 48 5 
lO ZD 4-o.J 


j_AC 7^ AO f 


f\ A CQ 

U.458 


c oo 
5.82 


oc m 
25.1)3 


FD 


LERG 


i n 


lOOl 1630+4534 


1 6 10 19 8 
iu ju jz.o 


4-45 14 96 


0.493 


6.06 


25.17 


J 


LERG 


1.26 


ItXJl 130/+3639 


1 1 07 97 07 


+16 19 16 4 

tju .jy iu.^t 


r\ CQQ 
U.583 


6.60 


oc in 
25.30 


J 


LERG 


1 1 A 

1 . 14 


7jh Alii i 1A 1 C 

/L 1)223+3415 


(1? 93 47 94 


+34 15 11 9 

TJ^ 1 J 11.7 


U.4/3 


c no 

5.y2 


oc cc 
25.55 


CD 


HERG 


/; ii 

6.31 


7P 17114_A£1S 


17 11 41 84 


4-66 18 56 7 

TUU do ju. / 


O 5A9 


f. 48 
0.45 


9C. A9 
2 j .02 


CD 


1 1 1 k 1 1 


j.o4 


7P 091 14-141 8 


09 1 1 98 1Q 


4-14 1 8 10 6 

TJt lo JU.U 


465 


5 87 


95 66 


V 1 > 


LERG 


fi 60 


TOOT 11014-1114 


1 1 01 1 9Q 


4-11 14 07 

TJJ J*+ U / .U 


565 

U.JDJ 


6 50 


95 66 


T 

J 


HERG 


1 06 


7r 091Q-I-1491 


02 19 37.83 


+34 23 11.2 


5Q5 


6 66 


95 Q8 


T 
J 


HERG^ 


6 30 

O.JU 


6P 08504-1747 


08 50 24.77 


+37 47 09.1 


407 


5 41 


96 1 5 

ZO. 1 J 


FD 


HERG 


Z.7J 


6r 1 9004-141 6 


12 00 53.34 


+34 16 47.3 


530 


6.29 


26.17 


CD 


LERG 


1.62 


6P 1 1 19-J-141Q 


11 32 45.74 


+34 39 36.2 


512 


6 1 8 


96 11 


^ i > 


HERG 


2 14 


6C 0857+3945 


08 57 43.56 


+39 45 29.0 


0.528 


6.28 


26.34 


CD 


HERG 


2.64 


3C 16 


00 37 45.39 


+13 20 09.6 


0.405 


5.41 


26.82 


CD 


HERG 


4.48 


3C 46 


01 35 28.47 


+37 54 05.7 


0.437 


5.66 


26.84 


CD 


HERG 


5.66 


3C 341 


16 28 04.04 


+27 41 39.3 


0.448 


5.74 


26.88 


CD 


HERG 


3.26 


3C 200 


08 27 25.38 


+29 18 45.5 


0.458 


5.82 


26.92 


CD 


LERG 


3.74 


3C 19 


00 40 55.01 


+33 10 07.3 


0.482 


5.98 


26.96 


CD 


LERG 


5.82 


3C 457 


23 12 07.57 


+18 45 41.4 


0.428 


5.59 


27.00 


CD 


HERG 


22.3 b 


3C 274.1 


12 35 26.64 


+21 20 34.7 


0.422 


5.55 


27.02 


CD 


HERG 


2.00 


3C 244.1 


10 33 33.97 


+58 14 35.8 


0.430 


5.61 


27.10 


CD 


HERG 


0.58 


3C 228 


09 50 10.79 


+14 20 00.9 


0.552 


6.42 


27.37 


CD 


HERG 


3.18 


3C 330 


16 09 35.01 


+65 56 37.7 


0.549 


6.41 


27.43 


CD 


HERG 


2.81 


3C 427.1 


21 04 07.07 


+76 33 10.8 


0.572 


6.54 


27.53 


CD 


LERG 


10.90 


3C 295 


14 11 20.65 


+52 12 09.0 


0.462 


5.85 


27.68 


CD 


LERG? 


1.32 



a Jet, Fat Double, Classical Double 
b Value obtained by Konar et al. (2009) 



TABLE 2 

Observation data for the ERA sample 



Source 


Instrument" 


Observation 


Exposure b 


Screened b 








time (ks) 


time (ks) 


TOOT 1301+3658 


C 


11568 


39.5 


39.5 


TOOT 1255+3556 


C 


11569 


39.5 


39.5 


TOOT 1626+4523 


C 


11570 


36.7 


36.4 


TOOT 1630+4534 


C 


11571 


20.8 


20.8 


TOOT 1307+3639 


c 


11572 


39.5 


39.5 


7C 0223+3415 


X 


551630101 


41.0/41.0/30.6 


32.5/35.7/21.4 


7C 1731+6638 


c 


11573 


24.0 


24.0 


7C 0213+3418 


X 


0551630201 


43.9/44.0/33.9 


15.4/27.1/3.4 


TOOT 1303+3334 


c 


11574 


39.7 


39.7 


7C 0219+3423 


c 


11575 


39.3 


39.3 


6C 0850+3747 


c 


11576 


39.2 


39.2 


6C 1200+3416 


X 


0551630301 


49.6/49.7/39.0 


37.6/41.0/24.6 


6C 1132+3439 


c 


11577 


39.6 


39.6 


6C 0857+3945 


X 


551630601 


24.9/24.8/18.2 


10.0/18.3/12.4 


3C 16 


c 


13879 


11.9 


11.9 


3C 46 


X 


0600450501 


17.9/17.9/13.2 


5.5/6.8/3.4 


3C 341 


X 


0600450601 


16.7/16.7/13.3 


13.3/16.1/11.2 


3C 200 


c 


838 


14.7 


14.7 


3C 19 


c 


13880 


11.9 


11.9 


3C 457 


X 


0502500101 


52.2/52.2/40.0 


36.8/39.2/28.8 


3C 274.1 


X 


0671640801 


27.2/27.1/22.6 


22.5/22.9/16.6 


3C 244.1 


c 


13882 


11.9 


11.8 


3C 228 


c 


2453/2095 


10.6/13.8 


24.4 


3C 330 


c 


2127 


44.2 


44.0 


3C 427.1 


c 


2194 


39.5 


39.5 


3C 295 


c 


578 


18.8 


18.2 



a Chandra, XMM-Newton 

b MOSl/MOS2/pn cameras for XMM-Newton 
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FIG . 1 . — On the left, 15 1-MHz radio luminosity vs redshift for narrow-line radio galaxies in the 3CRR, 6CE and 7CRS samples, with galaxies from the TOOT 
survey added in the range 0.4 < z < 0.6 (McLure et al. 2004). The McLure et al. sample is in the 0.4 < z < 0.6 redshift interval. On the right, the ERA 
subsample used in this project. 




FIG. 2.— 



The ERA sample; Radio luminosity coverage (left) and redshift coverage (right). 
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TABLE 3 

Details of radio maps used in the analysis 



Source Map frequency (GHz) Resolution(//) Ref. 



1 

1 

2 
2 
2 
2 
2 
2 
1 
2 
2 
1 
1 
1 
3 
3 
3 
3 
3 
3 
1 

4 
4 
5 
6 
3 



NOTE. — References: (1) Becker et al. ( 1995), (2) Mitchell (2005), 
(3) |http : / /www jb . man ac ■ uk/at las| (4) Mullin et al. (2008), (5) 
Hardcastle et al. (2002b), (6) Croston et al. (2005b) 



annuli around the source centroid. We used an annulus out- 
side the maximum detection radius to obtain the background 
counts, and these were subtracted from the source counts to 
obtain the net counts in each annulus. The point sources 
and radio emission identified during data preparation were 
removed and the annulus areas adjusted to account for the 
excluded regions. We used an energy range of 0.4-7.0 keV, 
this being within the reliable operating range for the Chan- 
dra data. For one source (3C 341), the extended emission 
was swamped by a bright nucleus, and we reduced the energy 
range to 0.4-2.0 keV to cut down the nuclear emission so that 
the profile could be modeled. 

For the XMM-Newton sources, since the pn camera is more 
sensitive than the MOS cameras, we created the pn profile 
first and then used the same annuli and background area for 
the MOS profiles. The three profiles were then scaled by their 
relative exposure times and added together. 

Table [4] contains the maximum detection radius and net 
counts within that radius for each of the sources. 

We fitted the surface brightness profiles with (3 models 
(see below) using the Markov-Chain Monte Carlo (MCMC) 
method described by Croston et al. (2008a) to explore the 
parameter space of these models and thus find Bayesian esti- 
mates of the core radii ( r c ) and f3 values. This MCMC method 
uses the Metropolis-Hastings algorithm in a manner similar 
to the METRO code by Hobson & Baldwin (2004), but im- 
plemented to run on a cluster of multiple processors using 
the Message Passing Interface (MPI: See Mullin & Hardcas- 
tle (2009) for more details). The method we use here differs 
from the the implementation of this method used by Croston 
et al. (2008a), in that we use a new fitting engine which al- 
lows the normalization of the (3 model and of any point-source 



component to vary freely during the fits. Plausible ranges for 
the free parameters were estimated and used to define uni- 
form (uninformative) priors for the MCMC method: for the 
normalizations and core radii (see below), priors uniform in 
log space were used to avoid bias towards large values. The 
uncertainties corresponding to 1 a errors for 2 interesting pa- 
rameters were determined using 1 -dimensional projections of 
the minimal n-dimensional volume that encloses 68% of the 
posterior probability distribution as returned by the MCMC 
code. This code was also used to determine the luminosities 
as discussed below. 

We initially used the appropriate instrument point spread 
function (PSF) alone to check whether this gave a satisfactory 
fit to the data. We then added a single f3 model to fit the ex- 
tended emission from the ICM (Cavaliere & Fusco-Femiano 
1976), convolved with the PSF. The surface brightness at ra- 
dius R is given by 

S(R) = S a (l + (£) 2 )-^+ - 5 
where r c is the core radius. Although the surface brightness 
distributions obtained from high-resolution / high sensitivity 
observations of many clusters have been found to be more 
complex than a single-fit /3-model, such a model provides an 
adequate parametrization of the surface brightness profile to 
integrate to obtain luminosity, given the comparatively poor 
data quality at these redshifts. The goodness-of-fit and f3 
model parameters are shown in Tableland Appendix lAlcon- 
tains the surface brightness profiles overlaid with the PSF and 
f3 model profiles. 

We could not obtain fits for six sources. TOOT 1630+4534 
could not be detected above the background (neither nuclear 
or extended emission); TOOT 1301+3658, 7C 0213+3418 
and 3C 16 had insufficient counts to create a profile, and 
7C 0223+3415 and 7C 1731+6638 had a point-source de- 
tection, but no extended emission above the level of the PSF 
wings. For these sources we estimated a 3a upper limit on the 
counts by obtaining the net counts within an estimated R500 
overdensity radius (see below). For the four faint sources, 
we used the median R500 of the 7C and TOOT sources; for 
3C 16 we used the median of the 6C and 3C sources. If the net 
counts were greater than three times the error on the counts in 
the background area, we used the net counts plus three times 
the error for the upper limit on the counts; otherwise we used 
three times the background error. (Since 7C 1731+6638 has a 
strong psf with no detectable extended emission, we used the 
background error method rather than the net counts method). 

The distributions of f3 and the core radius, r c , are shown 
in Figure [3] The majority of the j3 values are around 0.5, 
which is expected for groups and poor clusters (eg Hels- 
don & Ponman 2000), but there are three very high values 
of f3 (TOOT 1255+3556, 3C 46 and 3C 274.1 - these also 
have the three highest core radii) and two very low values 
(TOOT 1626+4523 and TOOT 1303+3334). The three TOOT 
sources are faint objects with low counts, so the model param- 
eters are very poorly constrained, and McLure et al. (2004) 
identified the host galaxy of 3C 46 as having undergone a ma- 
jor merger so its ICM may not be in hydrostatic equilibrium. 
3C 274.1, however, has an undisturbed host galaxy and so its 
steep profile is perhaps surprising. However, the profile fits 
the data well and since we are using the profile simply to ob- 
tain luminosity we are concerned only with the shape of the 
profile and not the accuracy of the (3 model. The j3 and r c 
values are degenerate and therefore not physically very mean- 
ingful, and the uncertainties on the extreme values are large 



TOOT 1301+3658 


1.1 


5.4 x 5.4 


TOOT 1255+3556 


1.4 


5.4 x 5.4 


TOOT 1626+4523 


1.4 


17 x 13 


TOOT 1630+4534 


1.4 


1.3 x 1.3 


TOOT 1307+3639 


1.4 


1.3 x 1.3 


7C 0223+3415 


1.4 


14.7 x 12.7 


7C 1731+6638 


1.4 


1.5 x 1.2 


7C 0213+3418 


1.4 


14.6 x 12.7 


TOOT 1303+3334 


1.1 


5.4 x 5.4 


7C 0219+3423 


1.4 


1.4 x 1.3 


6C 0850+3747 


1.4 


1.4 x 1.3 


6C 1200+3416 


1.1 


5.4 x 5.4 


6C 1132+3439 


1.4 


5.4 x 5.4 


6C 0857+3945 


1.4 


5.4 x 5.4 


3C16 


1.4 


1.2 x 1.2 


3C46 


1.5 


4.2 X 4.2 


3C341 


1.4 


1.3 x 1.3 


3C200 


4.9 


0.33 x 0.33 


3C19 


1.5 


0.15 x 0.15 


3C457 


1.4 


5.1 x 5.1 


3C 274.1 


1.4 


5.4 x 5.4 


3C 244.1 


8.4 


0.75 x 0.75 


3C 228 


8.4 


1.2 x 1.2 


3C330 


1.5 


1.5 x 1.5 


3C 427.1 


1.5 


1.8 x 1.1 


3C295 


8.7 


0.2 x 0.2 



AGN radio luminosity and cluster environment 



7 



TABLE 4 
Radial profile modeling results 



Source 


Rdet 


Counts a 


PSFfit 


PSF+/3 fit 





b 


r b 




kpc 




X 2 /dof 


X 2 /dof 






ai'csec 



TOOT 1301+3658 




(< 94) 




Low counts 


TOOT 1255+3556 


174 


340 


5.42/5 


0.23/5 


TOOT 1626+4523 


286 


108 


15.9/6 


0.617/6 


TOOT 1630+4534 




(< 63) 




Low counts 


TOOT 1307+3639 


162 


26 


4.18/4 


0.007/4 


7C 0223+3415 




« 77) 




Low counts 


7C 1731+6638 




(< 63) 




Low counts 


7L Oil J+J418 




(< 100) 




Low counts 


TOOT 1 '?n'?j.'?'^'^d. 


oo. y 


oO / 


^9 8/7 
3L.al I 


fi 1 zL/7 
0. 1+/ / 


/L, UZ1V+J4-ZJ 


1 f,A 
1 


Af, 


O.U/+ 


U.UUO/+ 


OL, Uojutj /^t / 




9^S 1 


o+.u/i' 


i n 9/q 


OL- IZULMOH-IO 




1 

JVoj 


1 ZU/ 1 u 


JO. 3/ 1U 




oov 


joy 


3 J.LI 1U 


q ^n/i n 

y.jy)/ 1U 


6L 085 1+3 747 


816 


612 


8.75/9 


2.17/9 


3C 16 




« 60) 




Low counts 


3C 46 


510 


170 


13.4/5 


0.57/5 


3C 341 


517 


214 


20.2/5 


0.25/5 


3C200 


200 


259 


16.1/6 


0.77/6 


3C 19 


589 


503 


294/8 


4.38/8 


3C 457 


1119 


2402 


51.4/9 


6.46/9 


3C 274.1 


610 


678 


121/8 


3.80/8 


3C 244.1 


966 


171 


27.8/6 


0.818/6 


3C 228 


537 


768 


32.5/9 


3.28/9 


3C 330 


473 


360 


82.1/10 


1.28/10 


3C 427.1 


675 


721 


288/10 


4.72/10 


3C 295 


719 


5308 


3978/17 


37.37/17 



0.49 
+0.02 
1.15 

30+ 19 
U ' JU -0.10 

0.49 

1.08 
0.16 
0.49 
0.49 
0.49 
22 +0.i3 
u - zz _0.05 
+0.54 
0.24 
+0.04 
0.02 
+ 0.03 
0.04 
n OO+0.07 

n 41 +0.15 

0.49 

+0.01 
0.97 
+0.96 
0.10 
+0.15 
0.05 
+0.06 
0.05 
+0.21 
0.11 
+0.34 
0.41 
+0.04 
0.05 
+ 1.10 
0.33 
+0.11 
0.07 
+0.14 



1.48" 1 



0.53] 



0.45] 
0.45^ 
0.41^ 



1.49] 
0.54j 
0.48^ 
0.61^ 
0.45^ 
1.16j 
0.4lj 
0.90j 
0.56 4 



5.58 

i9.oo±&° 

2.10±™ 
5.13 

1 qi+7.89 

5.24 
4.79 
5.29 

7.82±2:m 

4 7n+ 5 ' 28 
,u -3.70 

1 01+ 1 ' 49 
+0.07 
4.65 
+2.26 
6.73 
+ 15.80 
0.12 
5.74 
(-18.0 
-45.5 
+40.36 
0.23 
t-1.78 
-0.09 
+ 1.91 

I. 67 
+29.6 

17.3 

II. 5 
13.4 

+2.13 
0.01 
15.6 



6.22 
7.73 
1.12 

66.3 
1.74 
1.10 
7.73 
18.3 
32.0 
1.03 



0-40-0.02 

-, + 0.005 
-0.004 



0.50 1 



12 -31 6 . 2 
4.80± 2 ;^ 

3.1411;" 

r, cn + 0.19 
z - oz -0.15 



a Upper limits were obtained within estimated R500 ■ Counts for XMM-Newton sources are for the pn camera only 
b Italics indicate median values used for sources with low counts. 



and so cover more typical values. 

Luminosity was calculated by integrating the f3 model pro- 
file to the -R500 overdensity radius, using counts to flux con- 
version factors generated from the apec model. We extrap- 
olated the f3 model to the -R500 radius calculated using the 
R — T relationship from Arnaud et al. (2005): 

#500 = M^)~ 1 £Wir)' 3 , where h 2 (z) = fi m (l + zf + 

We calculated a luminosity for each sample of the output of 
the MCMC code, which provided us with a posterior proba- 
bility distribution function over luminosity, marginalized over 
all other parameters. We used the median rather than the mean 
of the posterior probability distribution function as our lumi- 
nosity estimate because the distributions were skewed for the 
fainter sources. Our quoted uncertainties on the luminosity 
are credible intervals defined on this one-dimensional pos- 
terior probability distribution function such that 68% of the 
probability is contained in the smallest luminosity range. The 
luminosity uncertainties take into account the (in some cases 
large) uncertainties on (3 and r c . 

Table [5] contains the X-ray luminosities for each source: 
within the 0.4 to 7.0 keV energy range to the maximum detec- 
tion radius; bolometric luminosity within the maximum detec- 
tion radius; bolometric luminosity within R500', and the bolo- 
metric luminosity scaled by h~ 1 (z) to correct for the critical 
density evolution. 

Nine sources have -R500 greater than the maximum detec- 



tion radius. For these, the extrapolated counts for the lumi- 
nosities within -R500 are less than 27% greater than the ob- 
served counts for all but one source - TOOT 1626+4523 - 
which has a shallow surface brightness profile and -R500 much 
larger than the maximum detection radius. Consequently the 
-R500 luminosity for that source is more than twice the lumi- 
nosity within the maximum detection radius. The mean ex- 
trapolation for the remaining eight sources is 12%. 

4.3. Spectral analysis 

Where possible, we used spectral analysis to obtain the 
ICM temperature, using the apec model for the thermal 
bremsstrahlung from the ICM and the wabs photo-electric ab- 
sorption model to take account of Galactic absorption. The 
high energy emission from the nucleus was excluded by re- 
ducing the energy range to 0.5 to 2.5 keV for all sources (0.5 
to 2.0 keV for 3C 341 - see Section l4~2l . In several of the 
XMM-Newton sources, a 1.5 keV instrumental aluminium flu- 
orescence line was visible, so we also excluded the 1 .4 to 1 .6 
keV energy band for all the XMM-Newton spectra (Freyberg 
et al. 2002). 

The spectrum was extracted from an annulus excluding 
the central nucleus and extending to the maximum detection 
radius, using the same background annulus as the surface 
brightness profile. We estimated the size of a region appro- 
priate for excluding the nucleus from the profile model fit by 
looking for the point at which the extended emission begins to 
dominate over the PSF Non-nuclear point sources and lobe- 
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FIG. 3. — Distribution of the f3 model parameters; /3 (left) and core radius (right). 



TABLE 5 
ICM X-RAY LUMINOSITY 



Source 


L x (0A-7.0keV) 


L x (bol) 


R500 


L x {bol) 


h(z)- 1 L x (bol) 




D ra d 


D ra d 


kpc 


R500 


R500 




xlO^erg.s- 1 


xlO^erg.s- 1 




xlO^erg.s- 1 


XlO 43 erp.s -1 



TOOT 1301+3658 






380 a 


< 1.26 


< 1.01 


TOOT 1255+3556 


25+ ' 10 


0.44t°;J° 


351 


0.48tg;? 7 


0.37 


TOOT 1626+4523 


1 61+ ' 70 

- L - D1 -o.99 


2.45+.^ 


533 


0,JU -4 .26 


4.15 


TOOT 1630+4534 






380 a 


< 2.44 


< 1.87 


TOOT 1307+3639 


23+ ' 11 
U - ZO -0.22 


40+ 20 


317 


41+ - 20 


0.30 


7C 0223+3415 




380 a 


< 1.45 


< 1.12 


7C 1731+6638 






380" 


< 2.97 


< 2.19 


7C 0213+3418 






380 a 


< 8.51 


< 6.60 


TOOT 1303+3334 


qp+1.89 
°- az -2.48 


fi 00+ 2 90 
°- uu -3.80 


525 




5.55 


7C 0219+3423 


56+ ' 24 
u - oo -0.53 


1 00+ ' 42 


380 


i.27t;-;s 


0.92 


6C 0850+3747 


qc+0.42 
z - ao — 0.49 


4 22+ ' 60 


648 




3.65 


6C 1200+3416 


8A8±\f 4 


12.69+i 55 


553 


a 9 o+0.95 
0-^0-0.92 


6.21 


6C 1132+3439 


a ci +1.81 

6.61I 200 


10.89+lp 


454 


' - JU -1.84 


5.54 


6C 0857+3945 


1.89li:S 


3-1^-2.35 


447 


O Q7+1-62 
z -°'-2.06 


2.16 


3C 16 






515 b 


< 2.79 


< 2.25 


3C 46 


3-90t?;^ 


5.95ti:i? 


535 


6.03+^-72 


4.79 


3C 341 


7O+0.34 


, 00+0.57 

i.oo_q S4 


429 


1 oq+0'51 
i - zy -0.54 


1.02 


3C 200 


1 41+0-42 


9a +0.69 
z - z8 -0.67 


474 


2 fi4+ - 93 
z - o4 -0.99 


2.07 


3C 19 


24.643+2 is 


36.589^323 


570 


36.4+ 2 f 


28.1 


3C 457 




9.03l|g3 


670 


D - 14 -1.13 


4.89 


3C 274.1 


„ 9r -+0.22 
z - zd -0.21 


or)-)-0.38 
o.oz_ 36 


343 


3.66t°;i 


2.92 


3C 244.1 


4.06+lH 


6.15±|;|| 


529 


4.541M 


3.61 


3C 228 


2.09t -j 8 


99 +0.67 
°- zz -0.89 


515 


99 +0.67 
°- zz -0.89 


2.39 


3C 330 


9 7 c-+0.42 


4 7K+0.72 
4 - '°-0.95 


430 


4 64+ - 70 


3.44 


3C 427.1 


I9.5lj;l 


27.9+31 


620 


26.2l 2 ; 5 


19.2 


3C 295 


125+J 


171±| ' 


966 


183t 4 


143 



a Median R 500 of 7C and TOOT sources. 
b Median R500 of 6C and 3C sources. 
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related emission were also excluded. We used a metallicity 
of 0.3 solar (Balestra et al. 2007) for all sources except for 
3C 295, 3C 330 and 3C 457, where metallicities had been 
calculated elsewhere (Belsole et al. 2007; Konar et al. 2009). 
These were 0.48, 0.2 and 0.35 respec tivel y. We discuss the 
effect of varying metallicity in Section IBTTl 

The ICM spectrum for the Chandra data could be fit imme- 
diately. For the XMM-Newton sources we first modeled the 
Galactic X-ray background in each camera using apec models 
for the two thermal components and a photon power-law with 
Galactic absorption for the extragalactic background. These 
were scaled to reflect the contribution from each camera, and 
the source spectra for the three cameras were then fit in the 
same manner as the Chandra sources, but with the three back- 
ground components. 

Once we had an initial estimate for the temperature, we var- 
ied the inner and outer radii to check that we had a stable value 
for the temperature with reasonable error limits, suggesting 
that the isothermal assumption was applicable across that an- 
nulus. For example, 3C 295 is known to have a cool core - 
Allen et al. (2001) found that the temperature rises with radius 
until about 10 arcsec. 

There were sufficient counts to model the spectrum for ten 
of the sources. For the fainter sources with a ft model fit, 
we extracted a spectrum and used the count rate from the j3 
model to estimate luminosity from the apec model. This was 
converted to temperature using the scaling relation of Pratt et 
al. (20 09), which provided an initial estimate of R500 (see 
Section l4~2l ). Iterating the process gave new estimates of the 
counts within R500, luminosity and temperature until the tem- 
perature converged. For the XMM-Newton sources, we used 
counts from the pn camera only in this process. The same 
process was used for the sources without a j3 model fit, using 
the 3cr upper limit on the counts. 

Table |6] contains the inner and outer radii of the annuli, the 
temperatures of the sources and the x 2 f° r the temperatures 
obtained by spectral analysis. As expected, the temperatures 
are for the most part typical of poor clusters with only three 
sources above 3 keV. The errors on most of the spectral tem- 
peratures are large, as is to be expected from such faint ob- 
jects. Temperatures have been obtained for some of the clus- 
ters by other researchers (Belsole et al. 2007; Allen et al. 
2001; Konar et al. 2009) and our values are compatible with 
their results. 

We have not included any corrections for the reported 
difference between temperatures obtained by Chandra and 
XMM-Newton since the difference at 3 keV and less is slight 
(Schellenberger et al. 2012). 

5. RESULTS AND DISCUSSION 
5.1. Radio galaxy environments 

In order to investigate the relationship between radio galax- 
ies and their cluster environments, we compared radio lumi- 
nosity and ICM luminosity for the full sample, for the differ- 
ent excitation classes and for the FR classes. These are plotted 
in Figures |4] and [5] 

The majority of the sources have similar ICM luminosity. 
Although the radio luminosity covers three decades, 80% of 
the sources have ICM luminosity within the single decade 
from 10 43 to 10 44 erg s -1 . The five sources outside this 
range are all LERGs (with the possible exception of 3C 295), 
spreading over about 2.5 decades, and X-ray luminosity does 
appear to rise with radio luminosity for the LERGS. Since we 



TABLE 6 
ICM TEMPERATURES 



Source 


Method 


Annulus radii 
kpc 


Temperature 
keV 


X 2 /dof 11 


TOOT 1301+3658 


Upper limt 




< 1.40 




TOOT 1255+3556 


Estimate 


15-174 


1 04+011 




TOOT 1626+4523 


Estimate 


15-296 


2 u +0.29 
z,± -0.81 




TOOT 1630+4534 


Upper limt 




< 1.69 




TOOT 1307+3639 


Estimate 


16-162 


_ „„_i_n 1 1 

0.98t{5;£ 




7C 0223+3415 


Upper limit 




< 1.41 




IV* 1 I J l+OO JO 


Upper limit 




*s 1 11 

< 1. II 




7C 0213+3418 


Upper limit 




< 2.37 




TOOT 1303+3334 


Estimate 


51^448 


q ■ 32 
62 




7C 0219+3423 


Estimate 


17-164 


1 36+° :is 
-0.91 




6C 0850+3747 


Spectrum 


21-160 


2 - 86 to'78 


3.77/4 


6C 1200+3416 


Spectrum 


252-944 


2 46+ 1 '' 89 


84.0/78 


6C 1132+3439 


Spectrum 


15-486 


71 +0.'69 
37 


4.16/5 


6C 0857 


Estimate 


251-816 


1 69+° :24 
i.oy_Q 55 




3C 16 


Upper limit 




< 1.78 




3C46 


Estimate 


170-510 


9 -, -, +0.14 
^ -0.20 




3C 341 


Estimate 


115-517 


1 '^°-0.22 




3C200 


Estimate 


15-215 


1 74+0.16 

1 - ' -0.23 




3C 19 


Spectrum 


29-383 


9 *y+0.63 

-0.38 


4.78/11 


3C457 


Spectrum 


168-839 


3 10+ 2 ' 95 


95.7/75 


3C 274.1 


Spectrum 


166-555 


oq5 +0.29 
U - MO -0 .25 


9.99/12 


3C 244.1 


Estimate 


14-966 


nc+0 .17 
z - uo -0.19 




3C 228 


Spectrum 


31-316 


O oo"l~2.39 
z,z -0.71 


1.06/4 


3C330 


Spectrum 


25^441 


1 61+ 1 ' 26 


0.76/3 


3C 427.1 


Spectrum 


32-193 


3 14+5-27 


2.86/4 


3C295 


Spectrum 


29-575 


6 09+ ' 73 


72.3/85 



a For temperatures from spectral analysis 



had censored data, we checked for correlations using general- 
ized Kendall's r tests (Isobe et al. 1986), and the results are 
shown in Table [7] We found a correlation between radio lu- 
minosity and ICM luminosity for the full sample and for the 
LERG subsample, but not for the HERGs. The LERG result 
does however depend heavily on 3C 295; the correlation is 
considerably weaker when 3C 295 is removed from the LERG 
subsample. 

We checked that this result was independent of black hole 
mass using the black hole masses calculated by McLure et al. 
(2004), who found a correlation between radio luminosity and 
black hole mass that was driven by the HERG subsample. We 
found no evidence of a correlation between ICM luminosity 
and black hole mass for the full sample or for the individual 
subsamples. 

Since the extrapolation of the luminosity to R500 depends 
on temperature, we checked the effect of propagating the tem- 
perature uncertainties through the calculations for some of the 
sources with large temperature errors and/or big differences 
between the maximum detection radius and R5Q0. For most 
of the sources, the difference was slight, but there were big 
changes to the uncertainties for some of the sources with large 
errors on the temperatures. 3C 427.1 had the biggest change, 
with the la range going from 28 — 30 x 10 43 erg s _1 to 
19 — 54 x 10 43 erg s _1 . This would however retain 3C427.1's 
position as one of the most luminous clusters and radio galax- 
ies and so the LERG correlation remains. 3C 457 has the 
next largest temperature errors and the upper limit on the lu- 
minosity potentially increases from 7 x 10 43 to 11 x 10 43 
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FIG. 5. — ICM X-ray luminosity coverage, separated into excitation classes. 

erg s _1 . Its position within Figure @] however, means that 
an increase in luminosity would merely weaken the HERG 
correlation further. Overall, we found that propagating the 
temperature uncertainties through the calculations would not 
affect the results of the correlation tests. 

A second potential source of error was our assumption of 
a metallicity of 0.3 solar for the majority of the sources. 
There were insufficient counts to allow metallicity to vary 
during modeling, so we recalculated the temperatures and lu- 
minosities of three of our sources with spectral temperatures 
(3C 427.1, 6C 1 132+3439 and 6C 1200+3416) with metallic- 
ities of 0.1 and 0.5 solar, these being outside the rms bounds 
given by Balestra et al. (2007). The biggest change in lumi- 
nosity was 4%, which was well within the luminosity errors 
of the sources. We also recalculated the luminosities for three 
of the sources with estimated temperatures. The biggest lu- 
minosity change was for the coolest source - TO 1255+3556 



- where dropping the metallicity to 0.1 solar increased the 
luminosity by 12%. Again, this was well within the la er- 
ror bounds. We therefore concluded that inaccuracies in our 
metallicities were unlikely to affect our results. 

As can be seen from the right-hand plot in Figure |4] any 
correlation for the FRII subgroup would be weak, and this is 
confirmed by the Kendall's r test. There is insufficient data to 
examine the FRI subgroup. 

5.2. Comparison with general cluster environments 

During the analysis, we used various assumptions based on 
the expectation that the cluster environments of our radio-loud 
AGN are not markedly different from other clusters of similar 
luminosity. In particular we assumed that the X-ray cluster 
luminosity is a good proxy for gravitational mass and that the 
ICM follows the expected luminosity-temperature relation. 
There is evidence that radio-loud groups of a given X-ray lu- 
minosity are hotter than similar radio-quiet groups (Croston 
et al. 2005a), but this effect is small enough not to be seen at 
higher luminosities (Belsole et al. 2007). 

We expected that our cluster environments, although rel- 
atively poor, would be sufficiently luminous for any disrup- 
tion by feedback to be smaller than the experimental errors. 
We checked the validity of this assumption by looking at the 
Lx — Tx relation, the gravitational masses and the entropy of 
the clusters. Figure |6]plots ICM luminosity (corrected for red- 
shift evolution) against temperature for the results obtained by 
spectral analysis, and Table|7]lists the results of the correlation 
tests. Since the majority of the temperatures were estimated 
from a scaling relation (Pratt et al. 2009), it would be startling 
if there were not a strong correlation for the full data set. The 
bulk of the scatter comes from the sources with temperatures 
derived from their spectra, and this is reflected in a weaker 
correlation for that subpopulation. 

Since our data are doubly censored, we used Schmitt's lin- 
ear regression (Isobe et al. 1986) to obtain the Lx — Tx re- 
lation. This is known to give a biased estimator when there is 
intrinsic scatter in both variables, so we took the bisector of 
the two regression lines (Hardcastle & Worrall 1999). 
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TABLE 7 

Generalized Kendall's t correlation tests 



Subpopulation Z 


N 


P 


Radio luminosity vs ICM luminosity 






All data 2.841 


26 


0.0045 


HERG 0.862 


15 


0.3887 


HERCj-no /C 0219 0.446 


14 


0.6557 


LERG 2.655 


11 


0.0079 


LERG-no3C295 2.140 


10 


0.0324 


FRII 1.903 


21 


0.0571 


Black hole mass vs ICM luminosity 






All J„,„ 1CTO 

All data 1.539 


26 


0.1237 


riLRLr 1.4/0 


15 


0.1415 


LERG 1.202 


11 


0.2292 


ICM temperature vs luminosity 






All data 5.082 


26 


< 0.00005 


Spectrum method 2.236 


9 


0.0253 


B 99 (977 kpc) vs ICM luminosity 






All data 1.965 


26 


0.0494 


HERG 0.963 


15 


0.3354 


LERG 1.456 


11 


0.1454 


FRII 2.064 


21 


0.0390 



B 9g (564 kpc) vs ICM luminosity 

All data 2.557 

HERG 1.878 

LERG 1.456 

FRII 3.064 



26 
15 
11 
21 



0.0106 
0.0604 
0.1454 
0.0022 



NOTE. — Z is the correlation statistic; N is sample size; p is probability 
under the null hypothesis 



The solid line in Figure [6] shows the regression line for all 
the ERA sources including the upper limts: Lx = 6.08 j^' 34 X 
10 44 (T X /5) 312±0 01 . The dashed line shows the Pratt et al. 
Lx — Tx scaling relation that we used to obtain the estimated 
temperatures, which is, as expected, very similar. 

Bohringer et al. (2012) predict a scaling relation slope of 
2.70±0.04 and their review of the literature cites slopes of 
2.6 to 3.7 for scaling relations derived from observations. Our 
result is therefore close to Bohringer et al.'s model predictions 
and compatible with existing observational results. 

We estimated the total gravitational mass for each cluster 
via the assumption of hydrostatic equilibrium, using equation 
(5.113) from Sarazin (1986). We obtained the density gra- 
dient from the (3 model parameters (Birkinshaw & Worrall 
1993). We then compared these results with the M—T scaling 
relation of Arnaud et al. (2005) and found them compatible 
within the lcr errors. 

We also calculated entropy S within O.I-R200 using 

h i ^ 3 (z)S = kT/n^ 3 , where -R200 is the radius at an over- 
density of 200 (Arnaud et al. 2005), kT is the ICM temper- 
ature and n e is the electron density. We converted (3 model 
counts to electron density via the volume emission measure 
from the apec model. We compared our results with the scal- 
ing relation of Pratt et al. (2010) and again found them within 
the expected range. 

We concluded that, within the experimental errors, the clus- 
ter environments of our radio galaxies are comparable to to 
those of other clusters of similar luminosity, and find no ev- 
idence (within our often large uncertainties, and bearing in 
mind that only ten of the cluster temperatures have been di- 
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FIG. 6. — ICM luminosity vs temperature for the temperatures obtained by 
spectral analysis. The heavy, solid line shows the Lx — Tx relationship from 
all the ERA results (including estimated temperatures and upper limits) and 
the dashed line shows Pratt et al. (2009)'s Lx — Tx relation. 

rectly measured) that the radio sources have altered the global 
cluster properties. 

5.3. Comparison with optical measures 

Herbert et al. (2013) have calculated the galaxy-quasar spa- 
tial covariance function B gq at two radii (977 kpc and 564 
kpc) for the McLure et al. (2004) sample, and we used these 
values to look for a scaling relation between ICM luminos- 
ity and Bgq (Figure [7] Table [7). There is only slight evi- 
dence of a relationship for the 977 kpc values, but the full 
564 kpc data set shows a correlation and this is strengthened 
when the FRI sources are excluded. Whether this stronger 
correlation is due to a difference in the two types of galaxy, a 
change in scaling relation with radio luminosity or some other 
cause cannot be determined from these data. For the FRII 
sources, Schmitt's linear regression gives a scaling relation of 
log lp (ix) = (0.0019 ± 0.0001) x B m + (43.28 ± 0.11), and 
this is plotted in Figure [7] 

Yee & Ellingson (2003) found a power law relation between 
B gc (the galaxy-cluster center correlation) calculated within 
500 kpc, and ICM luminosity - since their sample is more 
luminous than the ERA sample, they have no negative values 
of B gc . In order to compare our results with those of Yee and 
Ellingson, we removed our negative values from the 564 kpc 
data set and found a good correlation and a scaling relation of 



L x = l-711o:69 x lO^Bgq) 2 - 61 ^- 01 . Figure! shows a log 
plot of the B gq data overlaid with the two scaling relations. 
Our line is a little steeper than that of Yee & Ellingson; their 
sample contains brighter sources over a larger redshift range 
than ours and their cluster richness measures are calculated 
within slightly different radii, so a difference in the regression 
line is not unexpected. 

5.4. Implications 

Overall, we have found a correlation, significant at the 
~ 99.5% confidence level, between radio luminosity and clus- 
ter richness. This result is similar to that of Herbert et al. 
(2013), working with optical overdensity parameters of the 
McLure et al. (2004) sample, who found a 98% significant 
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FIG. 7. — ICM luminosity vs the galaxy-quasar spatial covariance function B gq (values taken from Herbert al. 2013). B gq is calculated within 564 kpc (left) 
and 977 kpc (right). The regression line is calculated for the FRII sources only. 




Bgq to 564 kpc 



FIG. 8. — ICM luminosity vs log 10 B gq for all positive values of B gq 
(calculated within 564 kpc) overlaid with the scaling relations from the ERA 
sample (solid line) and Yee & Ellingson (2003) (dashed line). Note that Yee 
& Ellingson used B gc calculated within 500 kpc. 

correlation between radio luminosity and environmental over- 
density within 564 kpc of the source. Thus for the sample as 
a whole, the optical and X-ray measures of environment rich- 
ness yield a similar result. 

Splitting the sample by spectral class, we have found a 
slightly reduced probability of correlation between radio lu- 
minosity and cluster richness for LERGs compared with the 
full data set, but no correlation for HERGs, suggesting that 
the overall correlation is driven by the LERGs alone. This 
adds to the body of evidence that there are fundamental dif- 
ferences in the properties of the two classes of radio galaxy. 
Our HERG clusters occupied a single decade of X-ray lumi- 
nosity (10 43 to 10 44 erg s" 1 ), while the LERG clusters had 
a much wider range of luminosities. However, there is about 
an order of magnitude of scatter in the relationship, and large 
uncertainties at the lower end of the luminosity range where 



most of the upper limits are also situated. This, taken with the 
small sample size, makes any conclusions tentative. 

Best (2004), Hardcastle (2004) and Gendre et al. (2013) aU 
found that HERGs lay within a limited range of relatively low 
environment richnesses. Hardcastle used B gg within 500 kpc 
to measure cluster richness. Bearing in mind the weakness 
of the correlation, the considerable scatter and the measure- 
ment radius of 564 kpc rather than 500 kpc, we can use the 
Lx — B gg relationship given in Section 15.31 to estimate that 
our clusters lie roughly within —150 < B gg < 400, which 
is roughly compatible with the range occupied by the HERGs 
in Hardcastle's sample. McLure & Dunlop (2001) and Har- 
vanek et al. (2001) report similar B gg ranges for their QSOs 
at z ~ 0.2 and z < 0.4 respectively. For their QSOs in the 
0.4 < z < 0.6 region, however, Harvanek et al. found higher 
values of B gg , roughly corresponding to cluster luminosities 
fromlO 43 to 10 45 ergs" 1 . 

Looking at results for higher redshift radio galaxies, Belsole 
et al. (2007) and Wold et al. (2000) both find HERG cluster 
luminosities within similar ranges. This could hint that any 
evolution of the HERG environment is slight, but it would be 
unwise to draw any conclusions based on this cursory analy- 
sis. We will look at evolution of the environment in the next 
phase of the program. 

Turning to the LERG subsample, the range of cluster lu- 
minosities is larger than for the HERGs, a result also re- 
ported by Hardcastle (2004), Best (2004) and Gendre et al. 
(2013). Best also found a 99.95% correlation between en- 
vironment richness and radio luminosity; our correlation is 
weaker (99.2%) but supports his result. This suggests a link 
between radio properties and the environment that does not 
appear to exist for HERGs. However, when 3C 295 is ex- 
cluded, the LERG sources are all within about 0.5 dex of the 
HERG range so the difference between the samples is now 
slight and the LERG correlation between radio and environ- 
ment luminosity becomes very weak. 

The answer to the question posed in the introduction: 
whether the radio-galaxy luminosity is related to large-scale 
cluster environment appears to be "Yes". There is, however, 
considerable scatter in the results, and the correlation is driven 
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by a small number of sources. We also found tentative ev- 
idence that it may be the population of low-excitation radio 
galaxies driving this relationship, with no correlation found 
for the high-excitation subsample alone. We have also shown 
that the correlation between Lr and environmental Lx (for 
the full sample) is not driven by a correlation between cluster 
X-ray luminosity and black hole mass. 

Such a correlation between radio luminosity and environ- 
mental richness would be expected if jet properties are de- 
termined by the properties of the surrounding hot gas; how- 
ever, there are many reasons why such a relationship may be 
expected to have considerable scatter (as we observe), even 
for the narrow redshift range we consider here. An important 
source of scatter is the still poorly understood relationship be- 
tween jet mechanical power and radio luminosity. Progress 
has been made in constraining this relationship observation- 
ally (e.g. Birzan et al. 2008, Cavagnolo et al. 2010, Godfrey 
& Shabala 2013), but there is one to two orders of magnitude 
scatter. This on its own could be enough to explain the scat- 
ter we observe in the relationship between radio luminosity 
and environment, even if jet power and environmental rich- 
ness were tightly correlated. It is worth bearing in mind, how- 
ever, that differences in environment are the most plausible 
origin for the scatter in the jet power/radio luminosity rela- 
tion, which would act in the opposite direction to tighten the 
relationship with environment. 

A second important source of scatter is the relationship be- 
tween central cooling time and total ICM X-ray luminosity. 
If jet power is controlled by the properties of the hot gas en- 
vironment, then it must be the central gas distribution that is 
important. Both cool core and unrelaxed clusters can be found 
at all luminosities, and the central gas density is uncorrelated 
with total ICM X-ray luminosity for the cluster population 
as a whole (e.g. Croston et al. 2008b). It has been found 
that low-power radio galaxies appear to require a cool core or 
dense galaxy corona (e.g. Hardcastle et al. 2001, 2002a; Sun 
2009); however, it has not been observationally established 
whether this is the case for the FRII population that form the 
majority of our sample. It is therefore plausible that there is 
substantial scatter between our measured X-ray luminosities 
and the central hot gas properties that may drive jet behaviour 
(at least in part of the sample), which could be a major con- 
tributor to the scatter we observe between Lr and Lx- 

6. CONCLUSIONS AND FURTHER WORK 

We have made a comparison of low frequency radio lumi- 
nosity and cluster environment richness of a sample of 26 
radio-loud AGN. We excluded any effects of environment 
evolution by taking the sample from a narrow redshift range at 
z ~ 0.5. The sample covered three decades of radio luminos- 
ity and contained both high- and low-excitation sources. Our 
measure of environment richness was ICM X-ray luminosity, 
obtained from Chandra and XMM-Newton observations. 

Our main findings are: 

• Over the full sample, there is a correlation between ra- 
dio luminosity and cluster richness, using ICM lumi- 
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nosity as the measure of cluster richness. There is how- 
ever total scatter of about one order of magnitude in en- 
vironment richness at a given radio luminosity, which 
is not much smaller than the total range of cluster rich- 
nesses. 

• There is tentative evidence for a difference between 
high- and low-excitation sources, with the HERGs oc- 
cupying a slightly narrower range of cluster richnesses 
than the LERGs and showing no sign of a correlation 
between radio luminosity and cluster richness whereas 
the LERGs have a similar strength correlation to the 
full sample. However, re-analysis without the bright- 
est source, whose classification as a LERG is disputed, 
reduces the correlation to only slightly above 95% sig- 
nificance. 

• We found no evidence for a correlation between radio 
luminosity and ICM luminosity for the FRII subpopu- 
lation. 

• Our results were compatible with published ICM 
luminosity-temperature scaling relationships. 

• We compared ICM luminosity and B gq in the hope of 
finding a scaling relation between the two cluster mea- 
sures. Although the large scatter makes any such rela- 
tion dubious, we found a correlation between log 10 Lx 
and B gq calculated to 564 kpc (Herbert et al. 2013) for 
FRII sources, and a power law relationship between Lx 
and B gq that is broadly compatible with that of Yee & 
Ellingson (2003). 

Having examined a sample of radio-loud AGN at a sin- 
gle epoch, we have found evidence for a correlation between 
radio luminosity and host cluster X-ray luminosity, as well 
as tentative evidence that this correlation may be driven by 
the subpopulation of low-excitation radio galaxies. This is in 
keeping with previous studies showing different accretion ef- 
ficiencies and host galaxy properties for the two types of radio 
galaxy. We also, in common with other researchers, found 
considerable scatter in the results, which may be a sign of 
more complex relationships between jet power and environ- 
ment than are generally assumed. 

During the next phase of the program we will compare our 
results with archive data and published results from different 
epochs to look for evidence of evolution of the environments 
of radio-loud AGN. 
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APPENDIX 

IMAGES AND SURFACE BRIGHTNESS PROFILES 

This appendix contains images and surface brightness profiles of the ERA sample in order of radio luminosity^. The images 
are of the X-ray emission overlaid with radio contours. The dashed circles are the maximum detected radius and the solid circles 
show R 500 . 

Twenty-two sources had sufficient counts to create surface brightness profiles. The PSF and f3 model profiles are overlaid on 
the profiles. Although profiles were generated for 7C 0223+3415 and 7C 0213+3418, they had insufficient extended emission to 
fit a P model. 



8 This version of the paper contains a subset of the images and profiles Astrophysical Journal article. These four examples show profiles with high 

in the full appendix, which will be available in the on-line version of the ./, . f , , v .,,, A , , . ,• 

rr and low counts tor Chandra and XMM-Newton observations. 
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FIG. 9.— 7C 0219+3423 




